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ABSTRACT 

Simple procedures are presented for treating cumulative fatigue 
damage under complex loading history using either the Damage Curve concept 
or the Double Linear Damage Rule. A single equation is provided for use 
with the Damage Curve approach; each loading event providing a fraction 
of damage until failure is presumed to occur when the damage sum becomes 
unity. For the Double Linear Damage Rule, analytical expressions are 
provided for determining the two phases of life. The procedure involves 
two steps, each similar to the conventional application of the commonly 
used Linear Damage Rule. When the sum of cycle ratios based on Phase I 
lives reaches unity. Phase I is presumed complete, and further loadings 
are summed as cycle ratios based on Phase II lives. When the Phase II sum 
reaches unity, failure is presumed to occur. No other physical properties 
or material constants than those normally used in a conventional Linear 
Damage Rule analysis are required for application of either of the two 
cumulative damage methods described. Illustrations and comparisons of 
both methods are discussed. 


INTRODUCTION 

Although the Linear Damage Rule (LDR) is commonly used in analyzing 
cumulative fatigue damage, it does not conform with the long established 
fact that order of loading significantly affects the summation of the cycle 
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ratios at failure. The non-conformity is particularly evident when high 
stress loadings are applied first, and low-stress loadings are subsequently 
applied until failure occurs. The sum of the observed cycle ratios is 
usually less than unity. Thus discrepancy is greater, the higher the 
ratio of the life at the low stress to the life at the high stress. On 
the other hand, if the low stress is applied first, followed by the high 
stress, the sum of the cycle ratios can be greater than unity. 

It was early in the history of the study of cumulative fatigue 
damage that the loading order effect was observed, and many approaches 
have been formulated to explain it. One of the explanations offered is the 
need for two linear damage rules. Cover [1] first suggested such a need on 
the basis of separating the fatigue process into one of crack "initiation" and 
"propagation", hypothesizing that the first linear damage rule related to 
initiation, the second to propagation. By further hypothesizing that the 
fraction of life required to initiate a crack was itself a function of 
life - being small at short life, and high at long life - it became easy to 
explain the order effect. Grover's treatment was qualitative, however, and 
he provided no quantitative formulation for separating the total life 
into its two components. 

Manson and co-workers followed Grover's lead in explaining the order 
effect [2,3], as well as developing an independent approach [4]. In [2] 
an attempt was made to quantify the "double linear damage rule" by providing 
an explicit formula for partitioning the total life into its initiation and 
propagation phases. Further study [3] led, however, to some question as 
to whether it is even valid to regard the two stages involved as "initiation" 
and "propagation". For a given loading the end of the "initiation" 
stage may be at one cycle ratio when this loading is followed by loading 
at one stress level, but could be at a different cycle ratio if followed 
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by loading at another stress level. In other words, "initiation" is a 
relative term, and depends on what is yet to come. Furthermore, careful 
examination of specimens loaded to a number of cycles, which, from previous 
tests, were known to be the effective "initiation" life, did 
not reveal the presence of measureable cracks. Thus, it could not be 
ascertained as to how large the crack must be to be regarded as "initiated". 
For this reason the terms "initiation" and "propagation" were abandoned and 
replaced by "Phase I" and "Phase II" to recognize that there was something 
different about the phases; although their physical interpretation was yet 
to be determined. 

In [3], Manson, Freche and Ensign studied the Double Linear Damage 
Rule (DLDR) behavior of a number of materials under two -block loading, and 
provided a formula for determining the "kneepoint" where the transition 
between the two phases occurred. Moderately good predictions could 
be made, but discrepancies were evident in some cases. It was concluded 
that foi" best practical results it was desirable actually to conduct the 
tests necessary to determine the kneepoint associated with the highest 
and lowest stress levels involved in a given application. The coordinates 
of the kneepoint then permitted the determination of two points on each 
curve of the two phases described by the DLDR, and it 
was suggested that the interpolated regions of the curves could be 
estimated by the investigator by sketching smooth curves through these 
points. No formulas were provided, however, thus making quantitative 
analysis difficult. 
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Since the publication of [3] attention has been directed by the present 
authors toward the solution of the two important remaining problems required 
for the practical implementation of the DLDR; an analytical procedure for 
determining the kneepoint without need for exreriment, and the development 
of equations for expressing the Phase I and Phase II life curves over the 
entire range of interest in a given application. It is the purpose of 
this report to explain the reasoning used to attack these problems and to 
present the solutions. 

RATIONALE OF APPROACH 

Although the final cumulative damage equations are simple, the 
process leading to their development involved several steps. The following 
discussion of these steps may be used in maintaining a perspective of 
why they were introduced. 

The first step is to establish a model for damage accumulation. It is 
recognized that the major manifestation of damage is crack growth. While 
the conventional treatment of fatigue involves tracking the growth of 
a single dominant crack, it is recognized that the early stages wherein 
such a dominant crack is developed involves many complicated processes 
such as dislocation agglomeration, subcell formation, multiple microscopic 
crack formation, and the independent growth of these cracks until they link 
and form the dominant crack. In order to account for these early processes, 
we make use of empirical information, established over a considerable period 
for many materials. Thus, we establish an "effective crack growth" 
equation which automatically accounts for the early processes without specifically 
defining or tracking them. This effective crack growth equation is then 
regarded as a model for damage accumulation. 
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Once the damage accumulation equation is established, either of two 
procedures can be followed. One is the direct use of the equation to track 
damage as loading is changed. We refer to this approach as the Damage 
Curve Concept (DC), which is already well-established from previous research 
(e.g., [5,6]). Whereas previous treatments used other functional forms 
of the damage curves, and parameters that required experimentation for 
each material, the method we have developed universalizes the parameters 
involved, and no experiments are necessary to implement the method. 

Because the DC concept results in rather complex numerical analysis, 
a second procedure is developed whereby the DLDR 
provides the basic framework. Rather than regarding the two 
linear regions as related to "crack initiation" and "crack propagation", 
they are regarded as two "phases" that result when the damage curves are 
replaced by two linear segments. Thus, the two phases are defined by 
their mathematical rather than their physical implications, although the 
two are loosely related, as will be discussed later. Replacing the DC 
procedure by the DLDR concept results in considerable simplification of the 
mathematics. 

The advantage of regarding the DLDR as a consequence of the linearization 
of the damage curves rather than resulting from the separation of the total 
life into its "initiation" and "propagation" phases is that it avoids 
the confusion associated with the observation made in [3] that "initiation" 
at a life level depended not only upon that life but also the later life 
level at which the fatigue process is continued. The result, as already noted, 
was the introduction of the terminology Phase I and Phase II to replace 
"initiation" and "propagation". However, in [3] it was not clear how 
to calculate the dividing point; hence an experimental approach was required. 
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A requirement for experimentation involves expense and is inconvenient 
when numerous exploratory calculations are required. Thus, a need has 
existed to determine the point of transition analytically. This goal 
has been accomplished in this report. 

By linearizing the generalized damage curves, the point of intersection 
of the two resulting straight lines is used as the transition. It will 
be shown that this point depends on the ratio of the lowest and highest 
life levels, N-j/N^, involved in the loading sequence to be analyzed. The 
analytical relation is determined not only from the linearized damage curves 
but from extensive experimental results obtained in previous investigations. 
Thus, the first step of an analysis is to determine from the loadings in 
a particular application the highest and lowest life levels involved, 
and then to determine from the formulas presented where the transition 
between Phase I and Phase II would be in a two-load-level test at the 
lowest life level and highest life level in the actual sequence being 
analyzed. That is, we determine how many cycles of the lowest life' level 
should be applied before changing over to the high life level so that all 
the Phase I life is consumed by the low-life loading, and all the Phase 
II life is consumed by the subsequent high-life loading. The concept 
involved here is analogous to the method of Corten and Dolan [7] wherein 
they replaced the regular S-N curve by a new one with coordinates 
established at the extremes of the life levels involved in the sequence 
to be analyzed. In the present method, however, we seek to determine 
two curves - one for Phase I and another for Phase II - in the range of 
the life levels entering the practical problem to be studied. Since from 
the transition point (kneepoint) we know the number of cycles to complete 
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Phase I for the lowest life level , we can get the life level for Phase 
II by subtraction from Likewise, since we know the number of cycles to 
complete Phase II for the highest life level N2, we can also get the Phase I 
life for N2 by subtraction. Thus, knowing the lives at the kneepoint 
provides us with two points on both the Phase I and Phase II life curves. 

If, indeed, the loading involved only two levels, knowing the 
kneepoint would be all that is needed to analyze any combination of these 
loadings. The discrete values of Phase I and Phase II lives at N-j and 
N2 are sufficient. If we assume that other loading levels intermediate 
between the two will also be present, then it is necessary to draw a 
curve for the Phase I and Phase II lives as a function of the loading 
parameter (stress, strain, load, etc.). We accomplish this purpose by 
choosing for the Phase I life an expression involving two parameters, 
so that knowledge of two coordinates will establish the parameters. The 
shape of the Phase I curve is also chosen so that if it is extrapolated 
to very long lives it will asymptotically approach the basic life curve 
(total life, consisting of both Phase I and Phase II) at very high 
life levels. This characteristic derives from the knowledge that for 
very long lives crack initiation occurs at a very high fraction of the 
failure life; thus we carry over the analogy between "initiation" and 
Phase I to display this characteristic as well. Using this restriction, 
a simple analytical expression is obtained for the Phase I life, and the Phase 
II life is obtained by subtraction from total life. 

Once the two life curves are known, the damage analysis proceeds 
just as with the LDR. However, it is performed in two steps. First 
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the cycle ratios are based on the life values from the Phase I curve. When 
such cycle ratios add to unity. Phase I is presumed to be complete. Then 
cycle ratios are based on the Phase II life curve. When the second sum 
reaches unity Phase II is complete and failure is assumed to occur. 

Although the procedure is developed by considering the coordinates 
of the kneepoint in a loading sequence with high load first, the resulting 
Phase I and Phase II curves are usable for any order of loading, even 
a two-level test with low load first. This can easily be seen in the 
specific examples shown later in the report. In fact, the procedure could 
have been devised by considering a low-high load sequence, using the 
correlations discussed in the report; the resulting Phase I and Phase II 
curves would be the same. Thus, these curves are usable for any 
sequence of loadings between the extremes. However, some small effects 
associated with loading order can be expected, as will be discussed in 
the report. 

In the following we shall develop both the DC and the DLDR concepts, 
and illustrate their application to some extensive experiments discussed 
in the literature. 

CRACK GROWTH MODEL 


In [8] we presented an equation for the number of cycles N qq^ required 
to develop a crack .003 in. deep in terms of the life to failure of a 1/4 in. 

dia. test specimen 


.003 


= Nf- 


2.5 N 


2/3 


( 1 ) 


In other reports [9,10] the number of cycles required to propagate a crack 
to failure after it had penetrated .013 in. was given as Thus, 
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the number of cycles to N required to develop a crack .013 in. deep is 

w = t\| _ 4N ^ (2) 

•^.ois ^ ' 

Since in both these cases the formulas were empirical, their forms being 
arbitrarily chosen to fit certain experimental results, it is possible to choose 
other analytical forms, keeping the numerical values by the altered 
formula to provide reasonable consistency with the numerical values implied 
by Eqs. (1 ) and (2) . 

For purposes of further application to the cumulative fatigue problem the 
form pursued was 

bn/ 

a = a^ + (a^ - a^)(ryN^) (3) 

In this formula N is the applied cycles to reach a crack length of a, and 
the characteristic defect length of the material when N^/N^ = 0. The 
value of a^ is the crack length at fracture when N^/N^ = 1. For a 1/4 dia. 
specimen it might be expected that a^ = 0.25, the full specimen depth. 

But, as pointed out in [2] the crack growth rate is so rapid near the fracture 
condition that almost any number close to 0.25 could be used. The value of 
a^ = 0.18 was chosen in [2] on the basis of several considerations. The 
equation then reduces to 

BN^“ 


a = a^ + (0.18 - aQ)(Ng/N^) 


(3a) 
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To obtain consistency with Eqs. (1) and (2) we first neglected a^, and 
recast Eq. (3a) by twice taking logarithms, resulting in 


Ln 


Ln(a/0.18) 

Ln(N^/N^) 


LnB + a Ln 


(4) 


Thus, choosing a = .003 and determining N from Eq. (1) for various values of 
N^, and then choosing a = .013 with N, determined from Eq. (2) we can plot 

T a 

the parameter on the left side of Eq. (4) vs. Ln N^, resulting in determination 
of B = 2/3 and a = 0.4. The crack length equation becomes 

a = a^ + (0.18 - (5) 

The numerical values predicted by Eq. (5) are quite close to those associated 
with eqs. (1) and (2) for the two crack lengths involved, assuming a^ = 0. 

For other values of a^, more appropriate values of a and B could be found, but 
for the present purpose this subject is not pursued further since, as 
will later be seen, the only number of significance in the cumulative 
damage application is a = 0.4. 


DAMAGE CURVE ANALYSIS 

Although this report is mainly concerned with the DLDR we divert 
briefly to a discussion of damage curves because their use aids us in 
determining the coordinates of the kneepoint of the cycle ratio curve 
needed for the DLDR. 

Cumulative fatigue damage analysis through the use of damage curves 
have been used by a number of investigators [5,6] as reviewed by Kaechele [11]. 
Figure 1 shows the basis of the approach. The curves plot the accumulation 
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of "damage" as a function of cycle ratio for various life values. The LDR 
would require that the curves be coincident for all life level, it is their 
separateness that produces the loading order effect. All curves start at 
the origin which represents the initial condition of the material where 
the damage state is zero, and terminate at failure F where D = 1.0. 

The DC concept is that damage accumulation proceeds along the curve 
associated with the life level at which a cycle ratio is applied. For 
example if a cycle ratio n.j/Ni = is first applied at the life level 

N-| , the damage will go from zero to A. If at this point a new loading 
level is introduced, the life of which is N2 and the damage curve for which 
is OBF we first find the point B at the same damage as A (since no 

damage is added by simply preparing to apply cycles at the next load level). 

If n 2 cycles are then applied at the N 2 life level, the point C is located 

as shown. Similarly, if n^/N^ is applied at the life level we locate 

the point D at the same level as C, and increment DE by the cycle ratio 
n^/N^, and so on for any arbitrary sequence of loading. When point F is 
reached failure occurs. 

The extent of the loading order effect is clearly shown in the figure. 
For example, in a two-level test at which the first load is applied at the 
low life level, say from 0 to A followed by a loading at a higher life level, 
say from B to F, the sum of the cycle ratios clearly omits the distance 
AB and is therefore less than unity. However, if the high life cycle 
ratio is applied first along OB, followed by the low life cycle ratio along 
AF, it is clear that the cycle ratio associated with AB is included twice, 
and therefore the summation of the cycle ratios is greater than unity. 
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Although the DC concept has been used by various investigators, there 
have not been any specific numerical assignments for the damage curve 
equations. Power laws for D such as (n/N^)“ have been used, where a is a 
function of N^, but specific functional forms for a in terms of the life 
values involved have not been proposed. For our analysis, we shall first 
derive equations based on the DC concept, using as a quantitative model 
for damage the analog of the crack growth Eq. (5). Thus, if we assume 

D= ^ a^ + (0.18 - aQ)(n/N^)(^/^^^f°‘'^ (6) 

we have a damage curve law that has the value D = 1 at n/N^ = 1 for all 
values of N^, and D = a^/0.18 is a small value for n/N^ = 0. 

Applying Eq. (6) in connection with the first increment of loading 
n^/N^ along OA in Fig. 1, we obtain the value of D at A and B 

Db = '’a = 08 [=0 ^ (“-'S - a„)(n,/N,)<2/3)Ni°-''] (7) 

Had we proceeded to point B along path OB, the cycle ratio required would be 
/^2 damage 
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Thus, equating Dg from Eq. (7) to Dg from Eq. (8), we get 

r\2'/\ = (9) 

When the added n 2 /N 2 cycle ratio is applied along BC, the total cycle 
ratio at point C is 

V^2 " + n2/N2 (10) 


If we proceed in the same manner to find the total cycle ratio at point E 
it becomes 


|(n^/Ni) 


(N^/N2)°-^ 


"2^^2/ 


(N2/N3) 


0.4 


"3/^3 


( 11 ) 


If point E coincided with point F, expression (11) would be set equal to 
unity, implying that failure occurs. In the more general case when K 
loadings are appl ied before failure occurs, the equation for DC analysis 
becomes 



Note that the subscripts 1, 2, 3, K-1 , K are the sequence numbers 

of the loadings as they occur, and in general will not coincide with 
increasing life level as used in the example of Fig. 1. While Eq. (12) 
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is very symmetrical, and its application is actually rather simple once 
the analyst becomes accustomed to its use, we will examine this method only 
cursorily in this report. Primarily, we shall make use of the specific 
equation for a two load level test in order to establish the kneepoint 
for use with the DLDR, as will be described later. Of particular interest 
in Eq. (12) is that the only constant that remains in the final equation 
is the exponent 0.4; all the others cancel in the derivation. In fact, 
even this constant is not critical. Examination of a series of other 
values in the range from 0.3 to 0.5 reveals the final results are not 
greatly altered compared to using the value 0.4. Thus, for the continued 
discussion, such use as will be made of Eq. (12) will be based on the 
universalized value of 0.4 as the exponent with the reasonable assumption 
that even if this value differs somewhat among materials, the results will 
not significantly be affected. 

SELECTION OF EXTREME LIFE LEVELS 

As noted earlier the implementation of the DLDR will be accomplished 
by selecting the constants involved in the calculations using the highest 
and lowest life levels of interest in a particular application. If these 
life levels are used to establish a kneepoint of a double-linear curve 
that results when the loading consists of only these two stress levels, 
two points will be established on the Phase I and Phase II curves. The 
complete curves are established by imposing additional reasonable requirements. 
All other loadings will then be analyzed as interpolations, using the same 
Phase I/Phase II curves. The question arises as to which are the two 
life levels to be so used. The lowest life level js not difficult to 
identify. It is usually associated with the grand hysteresis loop of the 
complete loading spectrum, and involves the total range from extreme-to- 
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extreme of loading. In many cases no mean stress will be present if 
sizeable plasticity is involved; if mean stress is present, however, its 
effect must be included to establish this lowest life, as will be discussed 
later. 

It would seem reasonable to assume that the lowest life level N, ,, used 

I ow 

to construct the kneepoint should have a certain minimum damaging effect, 
say 0.1% when calculated on the basis of a preliminary analysis using the 
LDR. For loadings leading to lives less than 1000 cycles this criterion 
would recognize the presence of even one cycle as significant. Of course, 
any high strainrange loading, if it influences the mean stress of later 
lower loadings, must be considered in this respect regardless of whether 
or not it is used to establish the kneepoint. 

The choice of the highest life of interest may present some 

difficulty. It should also be a life level for which some significant 
damage is present. For example, if a single cycle of loading at a 
condition that has a 10^ life is imposed, while all other loadings involve 

5 

lives less than 10 cycles, it will be clear in the final analysis that 
this single cycle imposes relatively little damage. It may then be 

5 

desirable to re-evaluate damage on the basis of the 10 life as the 
highest life involved, giving emphasis to those loadings that were most 
significant in the actual problem analyzed. Experience is a valuable 
tool. However, it should be emphasized that even though some loadings may 
be neglected in establishing the kneepoint, their damaging effects may 
still be included in the analysis, by extrapolating the resulting Phase I 
and Phase II curves beyond the extremes associated with the life levels 
chosen for establishing the knee. Some further discussion of this 
subject will be provided later but study is needed to establish criteria 
for extreme loading conditions. 
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ANALYTICAL REPLACEMENT OF DAMAGE INTERACTION BY ' .0 STRAIGHT LINES 
When only two levels of loading are involved Eq. (12) becomes 

(n^/N^) ' ^ + n2/N2 = 1 , or n2/N2 = 1 -(n^/N^) ‘ (13) 

If a plot is made of 02/112 vs. n.j/Ni a smooth continuous curve results 
passing through (0,1) and (1,0) at the extremities, as seen in Fig. 2. 

The DC concept is not based on the requirement of the representation 
of the results by two straight lines, and it is clear from the figure 
that no such straight lines result; the curves are continuous and 
of continuous slope. If, however, we wish to use these results as a 
basis for establishing an "effective" linear damage rule approximation, 
it would be necessary to replace the continuous curve by two straight 
lines. To obtain the coordinates of the kneepoint, i.e. the intersection 
between these two straight lines, we made use of several observations and 
required criteria: 

a) Since the portion of Eq. (13) involving N-j and N2 contains 
these terms within the ratio N-|/N2, then the coordinates of any 
effective kneepoint should also contain only the single parameter 
M-|/N2. Herein lies the major difference between the present method 
and that described in [2]. If the process involve’ is considered as 
"initiation" and "propagation", then the n.|/Ni coordinate should 
depend only on N-| , since it is at n-|/Ni that the crack "initiates". 

Once initiated, the crack is present for any level of N2 which is 
later used to propagate the crack. Similarly the vertical coordinate 
of the kneepoint should depend only on N2, since the cycle ratio 
to propagate the crack depends on the load level of the propagation 
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period, and not the previous load level to initiate it. By choosing 
both coordinates as functions of N^/N 2 we abandon the "initiation'V'pro- 
pagation" mechanism, and replace it with a Phase I/Phase II description 
comprising the linearization of the damage curves of Fig. 2. 

b) The choice of kneepoint should be made with consideration for 
the numerous two-level tests that have already been conducted and 
interpreted according to the DLDR. 

c) The single load level test must degenerate properly to the 
predictions of the two-load level tests when both load levels are 
equal. That is, in Fig. 2 the kneepoint must lie on the line joining 
(0,1) to (1,0) for N^/N 2 = 1. 

By first plotting the coordinates (n.|/N-|) and (n 2 /N 2 ) of the 
kneepoint for two-load level experiments versus N-|/N 2 on logarithmic 
coordinates we observed that resulting approximate straight Tines 
could be drawn having the same slope. Thus, we can represent the 
kneepoint coordinates by 


ni/Ni 


A(ryN2)^ 


knee 


Therefore 


02/^2 


= B(N^/N2)^ 


knee 


(14) 


[n^/N^ + r^^m2\nee = B)(N^/N2)^ (15) 

But when = N 2 = N, we know according to condition (c) above that 
ni + n 2 = N so that 


[(n-j+n 2 )/Nj = (N/N) = (A +B)(N/N)^ resulting in A + B = 1 (16) 
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A log-log plot of [n^/N^ + ^ 2 /^ 2 } 

knsG 

N 2 / N^ should therefore result in a straight line through the 
origin yielding the slope B. Such a plot is shown in Fig. 3, from 
which B = 0.25. 

Fig. 4 shows the plot of log (n«/No) vs. log (N,/N,) which 

^ knee ‘ 

according to Eq. (14) should yield a straight line. Constructing 
this line to be of slope 0.25 we obtain a value of intercept B = 0.65. 

Thus A = 0.35. The kneepoint coordinates are 

n,/N^] = 0.35(ryN2)°-^^ ; n 2 /N 2 ] = 0.65(N^/N2)°-^^ (17) 

knee knee 

Eq. (17) applies only when the life, , at the first loading level 
is shorter than the life, N 2 , at the second loading level and represents 
the relations of interest in establishing the Phase I/Phase II life 
curves. Implicitly, from these curves, the opposite loading order - 

long life loading followed by short life loading - automatically 
produces the appropriate result that the sum of the cycle ratios is 
greater than unity when only two loading levels are applied. However, 
explicit relations can also be obtained for this case by noting the 
symmetry of the damage curves in Fig. 2 about the diagonal when 
N 2 < N.J . The equations for this case can therefore easily be shown 
to be 

= 1-0.65 (ryN2)'°‘^^ ; =1-0.35(N^/N2) 

upper knee upper knee 


ni/Nij 


0.25 


( 18 ) 
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The effectiveness of Eqs. (17) and (18) in representing the 
kneepoints for tv/o-load level experiments conducted on four different 
alloys is shown in Figs. 5 to 8. The results for maraging 300 CVM 
steel and SAE 4130 steel shown respectively in Figs. 5 and 6 are 
from [3] and previously unpublished NASA data. Results for 

Ti-6A1-4V from [12] are shown in Fig. 7 and for D.T.D. 683 
aluminum [13] in Fig. 8. In each case the straight lines are drawn to 
the kneepoint determined from Eqs. (17) or (18) for the life levels 
used in the tests and shown in the figures. The dotted curves are the 
DC representation of the test conditions, using Eq. (13). In 
Figs. 5-7 the plots are for n 2 /N 2 vs. n-|/N^, while for Fig. 8 the 
sum of n.|/N.| and n 2 /N 2 is plotted against n.|/Ni since the data of 
[13] were presented in this form. Also of interest in Fig. 7 
are several plots for Ti-6A1-4V wherein a mean stress was maintained at 
the N 2 life level. The plots show that the kneepoint formulas 
are equally effective for this case as well as for completely 
reversed loading, although more data are required to verify this 
conclusion more generally. 

EQUATIONS FOR THE TWO PHASES IN THE DLDR APPROACH 
The importance of Eqs. (17) and (18) does not lie in analyzing 
two-load level tests; such loadings do not occur very often in service. 
The use is primarily for determining the life curves for the two 
phases in the DLDR approach to damage accumulation for more complex 
loadings. 
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The basic procedure for the DLDR is as follows: We assume that 

if the individual loadings within the block have lives from 
to = N2» any combination of such loadings will follow the same 

DLDR, with the kneepoint being determined from Eq. ( 17 ). The horizontal 
coordinate provides the value of n^ , where n-j is the Phase I 

life at the N-j life level. The coordinate n2/N2 provides 02 which is 
the Phase II life at the N2 life level, from which the Phase I 
life at N2 becomes (N2 - 1I2). Knowing two points on the Nj curve, 
we can then proceed by choosing the form 


Nj = exp (Z 


( 19 ) 


where Z and ij) are constants. This form has the advantage of containing 
only two adjustable constants, which can be determined from knowledge 
of the coordinates of the kneepoint. It also has the feature that 
if (j) is negative, Nj approaches as becomes very large. Thus 
the Nj curve does not intersect the Njr curve, but approaches it 
asymptotically as N^ approaches infinity, which correlates with 
our knowledge that for very long lives crack initiation occurs at a 
very large fraction of total life (if we accept an analogy between 
Phase I and crack initiation, but not their equivalence). 

To obtain two coordinates on the Nj curve we apply the equation 
for the kneepoint, Eq. ( 17 ) 


'I.N, 


= N, 


nj/Ni 


= 0.35Nj(Nj/N2) 


0.25 


■* knee 


( 20 ) 
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^I,N2 = ^2 


1 - 


= N, 


knee 


1 - 0.65(N^/N2) 


0.25 


( 21 ) 


Substituting into Eq. (19) allows the solution for Z and (t as 
follows 


Ln 


Ln 




Ln |o.35(N^/N2 


0.25 


Ln |i-0.65(N^/N2)°*^^ | 


( 22 ) 


Ln 


Z = 


0.35(N^/N2) 


0.25 


(23) 


Thus Eq. (19) is the relation for Nj, using the values of 4> and Z 
from Eqs. (22) and (23). 

The equation for Njj becomes 

Nji = Nf - Nj = [1 - exp(ZN/)] (24) 

As an illustration we consider a practical problem recently studied 
[12] involving a Ti-6A1-4V turbojet compressor disk at 77C(170F) wherein 
the loading cycle produced values of N-j = 2500 cycles and N 2 = 64,000 cycles. 
Using Eqs. (22) and (23), we obtain 4 > = - .523 and Z = -111.2, 
resulting in 

Nj = exp (-111.2 (25) 

Njj = N^[l-exp (-111.2 


(26) 
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Using the life curve as discussed in Appendix A, the results for 
this alloy are shown in Fig. 9. Appendix A, incidentally, illustrates 
an inversion procedure [14] which conveniently provides a closed 
form solution for when the basic life relation provided consists 
of an expression for strainrange in terms of life. 

APPLICATION OF THE LIFE CURVES 

Once the two life curves have been constructed, their use is 
analogous to the conventional linear damage rule. Cycle ratios are based 
on the Phase I life curve until the summation is unity; then cycle 
ratios are based on the Phase II life curve until again the sum is 
unity, at which time the region which is being analyzed is considered 
failed. If the life values used in the analysis are based on complete 
fracture of conventional (1/4 inch diameter) specimens in uniaxial loading, 
the type of failure implied in the analysis is the development of a 
large crack - of the order of 1/8 to 1/4 in. If baseline life values 
involve smaller size cracks (obtained either experimentally, or 
calculated as approximations using Eq. (5)) then failure can be 
defined as the development of other sizes of crack. In the 
following discussion we shall use only life data involving complete 
failure, leaving other refinements to later study. As an example 
of the application of the method, we refer to a recent study [12] of 
the Ti-6A1-4V compressor disk of a small gas-turbine engine. This disk 
contained a bolt hole which was analyzed as shown in Table I. A 
series of events for which strainranges and mean stresses were analytically 
determined, are listed in columns (2) and (3). The number of cycles 
of each event which occurred per loading block (mission) is listed 
in Column (4). Column (5) lists the cycles to failure for each 
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condition. Here the spline point inversion formula (Appendix A) was 
used, taking account of mean stress as outlined in Appendix B. Once 
is known for each condition the values of Nj and can be 
determined from Eqs. (25) and (26). These values are listed in 
Columns (6) and (7). 

To continue with the analysis is very simple. We sum cycle 
ratios based on lives of the Phase I process, as listed in Column (8). 
Since the sum of the cycle ratios based on Phase I loading is .0126 
per block, the number of blocks required to complete Phase I is 
1/.0126 = 79 blocks. Thus, after the first 79 blocks, damage is 
calculated according to cycle ratios based on Phase II loading, as 
shown in Column 9. Since the summation of cycle ratios based on 
Phase II loading is .005 per block, it requires 1/.005 = 200 blocks 
to complete this phase. Thus the total number of blocks is 70 + 200 = 
270 blocks. 

This calculation was checked by subjecting specimens to the block 
loading history. In 220 blocks, cracks of about .015 to .030 inch 
developed. Since the calculations are based on complete fracture 
of 1/4 in. dia. specimens, they naturally predict "failure" after 
a larger number of blocks, so the agreement must be regarded as very 
good (We are currently studying the question of modifying the value 
of to use life to a particular size of crack, rather than complete 
rupture of test specimens, but since such data are more scarce than 
fracture data, our decisions are not yet firm). By comparision, the 
LDR for this problem predicts 306 blocks, which also is reasonably 
accurate in this case. The reason for this result is that there 
are relatively few cycles in each block of the loadings that produce 
small damage per cycle; thus these loadings produce a relatively 
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small amount of the total damage. Increasing their damaging effect 
by interaction with the high amplitude loadings does not have the 
dramatic reduction effect on life that it does at the knee of a 
two-level test when the actual number of N2 cycles is very large 
relative to the N.j cycles. In other cases, - for example, during 
low amplitude vibration following high amplitude loading, - the 
life reduction predicted by the DLDR would be much greater compared 
to the LDR. 

DISCUSSION 

That Eq. (17) is satisfactory for determining the coordinates of 
the kneepoint is amply demonstrated by Figs. 5-8. Other tests described 
in the references from which these figures were drawn provide additional 
verification, but are omitted here for brevity. Thus, it can be 
concluded that these formulas are adequate for representing 
two-level fatigue tests. However, the usefulness of the approach 
must be judged on the basis of whether more generalized loading sequences 
can be properly analyzed. Unfortunately the test programs for which 
data are available were not designed to test critically the difference 
between the two methods. In the program results to be discussed, the 
computed lives by the two methods differed only moderately. Nevertheless, 
they provide a useful framework for discussion. 

The Three-Level Test Program of Webber & Levy 

In [13] Webber and Levy describe a very extensive fatigue 
test program on an aluminum alloy. Their main interest was to study 
data scatter in cumulative fatigue, so many tests were conducted at 
three stress levels, both singly and in sequences involving two 
stress levels, or all three. Statistical analysis of their data 
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provided information on fatigue at any chosen level of probability of 
failure. For brevity only baseline fatigue properties at each 
stress level, and experimental results for multiple-stress tests, 
will be shown for the 50% failure probability of many duplicate 
tests, thereby permitting data scatter to be ignored. 

The main results of the analysis are summarized in Table II. As 
noted, the three life levels involved were 51,900, 414,140 and 13,800,000 
cycles. The two-stress level tests have already been summarized in 
Fig. 8. The results shown in Table II refer only to three-level tests, 
each duplicated 15 times. Each of the nine three-level tests consisted of 
repetitive blocks of 9500 cycles within which loadings occurred at 
all three stress levels. The percentages of the loadings at each 
stress level are shown in Columns 2, 3, and 4. Thus, for example. 

Test 1 consisted of 95 cycles (1% of 9500) applied at the stress which 
alone yielded a life of 51,900 cycles, 5985 cycles at the stress for 
414,1400 life level and 3420 cycles at the stress which lasted 
13,800,000 cycles when applied alone. The observed average (assuming 
a log-normal distribution) cyclic life for each set of loadings is 
listed in Column 5. Analysis by the LOR resulted in the cyclic life 
shown in Column 6. For the DLDR two calculations were made, one 
using N-| = 51 ,900 and N 2 = 13,800,000 for the determinations of the 
kneepoint and the other using N-j = 51 ,900 and N 2 = 414,140. 

The damage of the cycles at the 13,800,000 life level was then 
included by extrapolating the Phase I/Phase II life curves beyond 
those used to establish the kneepoint. The results for both calculations 
are shown in Columns 7 and 8. Damage Curve analyses were also made 
according to Eq. (12), using the two exponents, 0.40 and 0.33, as shown 
in Columns 9 and 10. Normalization of the calculations of Column 
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6, 7, and 9 relative to the experimentally determined lives are 
shown in Columns 11 to 13. 

Features of the Methods as Revealed by Analysis of 
Experimental Programs 

We can now use results of Tables I and II as a framework around 
which to conment on certain features of the methods. 

Treatment of mean stress . In evaluating the currently formulated 
DLDR we note that the input required is identical to that ordinarily 
required for any cumulative damage analysis; no additional data what- 
ever are required . By expressing Nj and Njj in terms of using 
Eqs. (19) and (24) we can account for mean stress as easily as 
not. When mean stress is present, is calculated according to any 
of several conventional methods which take account of the mean 
stress, and the Nj and values for this condition involving 
mean stress are determined as the intersections of the Nj and Njj 
curves along the horizontal passing through the associated value of 
N^. One way of accounting for mean stress is discussed in [12], 
and illustrated in Appendix B. 

The manner in which the transition is made from one stress level 
to another can be very important in relation to mean stresses developed. 
Minor nuances of detail can lead to different mean stresses, particularly 
when load levels in the high cycle fatigue range are mixed with low 
cycle fatigue loading. This subject is discussed in [12]. But it 
is important to recognize that an appropriate hysteresis loop analysis 
is required before a cumulative fatigue damage analysis can be under- 
taken. 

Comparision of Double Linear Damage Rule with Linear Damage Rule . As 
can be seen by comparing Columns 5, 6, and 7 of Table II the DLDR 
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generally produces better life predictions than the LDR. On 
average the DLDR estimates life for this service approximately 125^ 
high (Column 12) while the LDR provides estimates 37% too high. In 
general, it can be expected that the DLDR will provide lower estimates 
of life because it allows for the damaging effects of high loadings 
on the life at subsequently applied lower loadings. Unusual cases 
involving only progressively increasing loading may, however, result 
in longer lives according to the DLDR than LDR, as suggested by Eq. (18) 
and illustrated in Fig. 5 by the test involving N.| = 584,750 and 
N 2 = 990 cycles. The degree of difference resulting from calculations 
by the two methods depends, however, on the relative cycle ratios 
involved in the loading. As can be seen in Figs. 5 and 6, large 
differences exist if the cycle ratios involved are near the kneepoint, 
in some cases by as much as a factor of 5 or more. However, if the 
conditions involved are .such that o.|/N.| or n 2 /N 2 are close to unity, 
the other cycle ratio is nearly the same by the two methods. In other 
words, if nearly all the damage is done by one of the loadings, the 
contribution of the other loadings is small regardless of which method 
of calculation is used. This, in fact, is illustrated by both the 
experimental programs of Tables I and II. For the test program of 
Table II, in each case the cycle fractions associated with the 13,800,000 
loading is very small, and the calculated damage contributed by this 
loading is relatively small. Thus, these tests are basically two-step 
loading tests for purposes of cumulative fatigue damage analysis. The 
combinations of loading are such that the operating condition of the 
414,140 cycle loading is dominant and it combines with an n.j/N.| loading 
which is remote from the kneepoint. 
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Thus, neither test program is designed to accentuate possible differences 
that can be calculated to occur according to the two methods. But 
sufficient differences can be seen in both cases to demonstrate that 
the DLDR is generally more accurate that the LDR. 

A significant implication of the DLDR relative to the LDR is 
that high loadings can drastically lower resistance to later lower loadings. 
This can explain why "endurance limits" can be so greatly lowered in 
service applications involving some high loadings. It also suggests 
that vibrations of low amplitude involving large numbers of cycles 
when superimposed on large amplitude loadings can be very damaging 
even though their presence might be adjudged innocuous according to 
the LDR. 

Comparision of Double Linear Damge Rule with Damage Curve 
approach . The DLDR approach is intended as a simplification of the Damage 
Curve concept, and predictions by the two methods will generally be 
reasonably consistent. This is shown in Figs. 5 to 8 for the two-level 
tests, and by a comparison of Columns 9 and 10 with Columns 7 and 8 in 
Table II. Whereas programming Eg. (.12) for Damage Curve analysis can 
be accomplished easily on a computer, or even hand calculator when 
only a few loadings are present, it can become cumbersome when 
numerous loadings, such as those of Table I, are involved. Since 
the results are usually quite similar by the two methods, the DLDR 
may be regarded as preferable. 

It should also be noted that whereas the Damage Curve method 
recognizes every change of loading order as significant in the 
calculation procedure (whether or not numerically significant in 
result), the DLDR does not place significance on loading order 
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within each of the Phases. Thus, once the transition from Phase I 

and Phase II has been established, rearranging the loadings within 

Phase I or Phase II in orders other than as they actually occur does not 

alter the result of the life calculations. For example, in Test 1 

of Table II DLDR calculations show that the transition from Phase I to 

Phase II occurs after 17.3 blocks of 9500 cycles, and Phase II is completed after 

the subsequent 36.2 blocks. Thus, in Phase I there are 1644 cycles 

of the 51,900 life level, 103,540 cycles of the 414,140 life, and 

59,166 cycles of the 13,800,000 life. In reality they occur in 

regular alternations. However, suppose we applied all the 1644 

cycles of the 51,900 cycle life first, then all the 103,540 cycles of 

the 414,000 life, followed by the 59,166 cycles of the 13,800,000 

life. According to the DLDR the life would be the same since all 

the loadings within Phase I are the same, and Phase I lives are the same, 

and a linear damage rule is applied within this loading period. Does 

it not, actually make a difference? As an indication of what to 

expect we made a calculation in which all the Phase I and Phase II 

loadings were applied in decreasing load sequence, and another 

calculation in which they were applied in increasing load sequence. 

By the DLDR the results were, of course, unaffected, but by the 
Damage Curve approach differences were obtained. 

In the decreasing load sequence, failure, according to the Damage 
Curve approach, would occur before all the loadings could be applied, 
and we would lose 30,970 cycles of the 414,140 loading in Phase II and 
123,894 of the 13,800,000 loading. Thus the loss is 30,970/41,140 + 
123,894/13,800,000 = 8.4% of the life fraction calculation by the 
DLDR. However, if the loadings are applied in decreasing life 
sequence, failure would not occur when the loadings allowed by the 
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DLDR had been imposed. Additional life would be possible, the magnitude 
depending on whch stress level was used to complete the test. If 
the 51,900 life level were used, the remaining life would be 11,700 
cycles, or 22.6% of this life level. If the 414,140 life level were 
used, the remaining life would be 43,720 or 10.5% of the life level. 
However, if the 13,800,000 life level were used the remaining 
life would be 373,500 cycles, or 2.7%. Thus, although the two methods 
indicate a difference in result, the effect is not very large. In 
service, interspersing of loadings is the more common occurrence, making 
the two methods agree more closely. But the subject should be 
studied further to gain insight as to what actually does happen, 
and to evaluate the relative accuracy and desireabil ity of the two 
methods. 

Choice of life levels for determination of kneepoint . The two life 
levels used to establish the kneepoint will generally be the highest 
and lowest life levels significantly involved in the loading history. 

In this way the Phase I and Phase II life curves are established over 
the complete loading history involved. However, ambiguities may 
develop as to whether the existence of relatively few cycles of high 
cycle loading qualifies that high life for use as the N 2 for the entire 
loading sequence. To gain some insight into this question the 
calculations shown in Column 8 of Table II were performed. In these 
calculations the = 13,800,000 cycles were ignored in establishing 
the kneepoint. N.| was taken as 51,900 and N 2 = 414,140 to establish 
the Phase I/Phase II life values. However, these curves were 
extrapolated to determine the Phase I and Phase II lives at 13,800,000 
cycles to failure, so that the damage due to this loading could be 
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included. Column 8 of Table II shows that in most cases the results 
differed little from those of Column 7 wherein the highest life of 
13,800,000 cycles was actually used to establish the kneepoint. For 
test 6, however, wherein the largest percent (81%) of the block 
was at the 13,800,000 life level, an 11% difference in life prediction 
occurred, but for such a loading sequence ambiguity would not normally 
arise in the choice of N 2 because the 13,800,000 cycle loading occurred 
for 3% of life during the loading history. Thus it can be concluded 
that when ambiguity occurs the high cycle loading can safely be 
ignored in establishing the kneepoint. But there is no reason to 
omit its damaging effect. By extrapolating' the Phase I and Phase II 
life curves to the long life required, the damage can be included. 

Since the two life curves are tied to the curve by the well- 
conditioned Eq. (19), there is little chance that extrapolations will 
diverge far regardless of any reasonable choice of N 2 in establishing 
the kneepoint. 

Physical significance of the two phases . One of the major 
differences between the DLDR approach described in this report and that 
proposed in [2] and [3], is the abandonment of the terms "initiation" 
and "propagation" as the physical description of the two phases. When 
these descriptors are used the implication has to be that n-|/Ni at the 
kneepoint must be a function only of N.j , and n 2 /N 2 at this point a 
function only of N 2 * Figure 10 shows that the values of n-i/N-j at the 
kneepoints for various tests in [2] and [4] do not correlate with 
N^. Furthermore, if the crack lengths at the kneepoint are calculated 
from Eq. (5) (letting a^ = 0), the values obtained for the tests 
conducted are so small as to be undetectable under normal circumstances. 
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The kneepoints for most of the tests were at values of n.|/N.j less 
than 0.2 or 0.3, and the calculations indicate that detectable cracks 
of the order of .003 in. require cycle ratios n-|/Ni greater than 
0.8 to 0.9, thus explaining why our early attempts to detect cracks 
at the kneepoint were unsuccessful. Thus, treating the DLDR as 
a process of initiation, followed by a propagation phase amenable 
to Linear Elastic Fracture Mechanics (LEFM) would involve applying 
LEFM to extremely small cracks, which is questionable. 

Our approach of determining the coordinates of the kneepoint as a 
function of N-|/N 2 not only produces better correlation with experiment, 
but may have physical significance in relation to the type of 
damage involved and why the perception of the end of Phase I at the 
N-j life depends on N 2 * For example, we can assume that the early 
stages of fatigue consist of the development of damage sites 
associated with local plastic deformation. In Phase I these sites 
are discrete and separated from each other. Phase II produces 
linking and the development of a dominant crack. Thus, how the two life 
level loadings interact depends on how the second loading contributes 
to the transition from one mode of damage accumulation to the other. 

For this reason both loading levels are involved in the transition. 

This view is consistent with that expressed in [7]. More study is 
obviously needed to clarify the physical picture, but it is fortunate 
that the engineering treatment presented does not require a detailed 
knowledge of the physics. The formulas provided in this report 
permit engineering calculations even though the full understanding of 
the physics lie ahead. 
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SUMMARY AND CONCLUDING REMARKS 

Although a rather detailed discussion has been provided to 
explain the rationale for the DLDR approach the practical procedure 
is extremely simple: 

a) Choose the lowest and highest life levels N^ and N 2 involved 
in the loading sequence of interest, the choice not being critical. 
From these two life levels, and using Eqs. (19) and (24), calculate 
the life curves for Phase I and Phase II. 

b) Analyze loadings just as would be done in conventional Linear 
Damage Rule analysis, except do it in two stages. Instead of 

using the N^ curve to establish cycle ratios, use the Nj curve and add 
cycle ratios in the conventional manner until the summation reaches 
unity; Phase I is now complete. Then for continued loadings, sum 
cycle ratios using the Njj life curves; when these ratios sum to unity, 
Phase II is complete, and the part is considered "failed". 

Thus, application of the DLDR requires no input information other 
than that required for a conventional Linear Damage analysis. Its 
main purpose is to account for a loading order effect known to occur 
with smooth specimens, and particularly to account for damaging 
effects of high loadings followed by low ones. The method does 
not pre-empt need for detailed analysis of hysteresis loops to 
establish strainranges and mean stresses for each event, which is 
also basically required when a LDR analysis is made. 

Although a number of practical cases have been studied by the 
DLDR, experience is still limited. In general it can be expected that 
when high and low amplitude loadings are interspersed, predictions 
made by this method will be more conservative and accurate than 
those made by a Linear Damage Rule, which is desirable to promote 
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safety. However, more study is needed to establish its true 
validity and value. More study is also needed to obtain a good 
perspective of its merits relative to the Damage Curve approach which 
in concept can be applied easily, although requiring laborious 
analysis for histories involving numerous events. 
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APPENDIX A - ILLUSTRATION OF THE LIFE INVERSION PROCEDURE 
The inversion procedure is described in detail in [14]. Basically 
the 1 ife relation is given in the form 

^ = e): (2N^)^ + ^ (2N^)*^ (A-1) 

To illustrate the procedure we refer to a recent study [12] on Ti-6A1-4V 
at 170F for which the equation is 

2.852 (2N^)‘-^°^^ + 0.01987 (A-2) 

In the spline point approach to the inversion procedure we first 
determine the transition strainrange and life, which in this case are 


AE.J. = = 0.0181 (A-3) 

where E=1 12x10^ MPa (16. 27x10^ ksi) 

Hj = [l/2 = 291 cycles (A-4) 

We then calculate three constants 

6 = (c/b)*^-^^ = 1.7704 (A-5) 

a = -0.17-0.52 Ln(c/b) = -1.2071 (A-6) 

p = 0.5(c/b)-0.3 = 2.2743 (A-7) 

The inverted life relations become 

N^ = Ny(Ae/Aej)^^^ exp i6(Ae/Aej)“j for Ae > ACj (A-8) 

= 3.43 exp(,01396AE'^-^°^b forAe>.0181 (A-9) 



36 


and 

= Nj(Ac/Ae-|-)^^^ exp6 (Ae/ACj)^ , for Ae < AGj (A-10) 

= 1.953 X exp(1.34xlO®Ae^*^^^^) for Ae<.0181 (A-ll) 

The agreement in life representation by Eqs. (A-9) and (A-ll) compared 
to representation by Eq. (A-2) is extremely close [12,14]. 
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APPENDIX B - MEAN STRESS EFFECTS 

The life relation such as Eq. (A-1) is valid for completely reversed 
strains without presence of mean stress. In the more general case when mean 
stress is present the approximate equation is in the form 

E -— (B-1) 

where k is an empirical constant [14] which should be determined experimentally. 
Since its introduction by Morrow (with k=l ) [15], Eq. (B-1) has been 
used with considerable success in the long-life nominally elastic region. 

It has in fact been extended to the plastic regime by a number of investigators, 
e.g. Landgraf [16] and Manson [14]. Under elastic conditions the plastic 
strainrange is negligible, and its relation to the elastic strainrange 
is of little practical significance. However, at larger strainranges, 
wherein the plastic strainrange is an appreciable part of the total, it is 
important that the relationship between the elastic and plastic strainrange 
be maintained so as to imply a reasonable cyclic stress-strain curve in 
the presence of mean stress. Since Eq. (B-1) implies large differences 
in the cyclic stress-strain curves, depending on whether tensile or 
compressive mean stress is present, it is desirable to reconsider the 
appropriateness of retaining the same plastic line in the presence of mean 
stress, which is implied by Eq. (B-1). 

Consider two cases: in the first tensile mean stress is present; 

in the second the mean stress is equal in magnitude, but compressive. 

If an alternating strainrange is applied in both cases, it would 
be expected that the associated stressranges would be equal, since the 
individual stresses in the tension case and compression case would be 
approximately equal but of opposite sign. In other words, the 
cyclic stress-strain curves should be ne&rly the same for tensile or 
compressive mean stress. One way of leaving the cyclic stress- 
strain curve unaltered by the presence of mean stress is to assume that 
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the transition strainrange is unaffected by mean stress. Then the 
elastic line would be shifted according to Morrow's model, and the plastic 
lines would also shift to retain the same transition strainrange. When the 
transition strainrange is the same, the ratio of elastic strainrange to 
plastic strainrange is dependent only on strainrange; thus the cyclic 
stress-strain curve is unchanged by presence of any mean stress. It is 
possible, however, that, in a second case, mean stress does alter the cyclic 
stress-strain relation. In order to retain symmetry between the curves 
for both tensile and compressive stresses, the assumption should be made, 
therefore, that the transition strainrange depends only on the absolute 
value of mean stress, not its sign. Then, the relation between elastic 
and plastic strainrange wiTl still depend only on strainrange, and the cyclic 
stress strain curves will be the same for equal values of both tension 


and compression, although different from the basic cyclic stress-strain 
curve for zero mean stress. 

In order to generalize the possible behavior, we can assume that 
in the presence of mean stress the transition strainrange changes from 
Aej (at zero mean stress) to k'Aej with mean stress, where k' depends 
only on absolute value of mean stress. Then, since the transition life 
must lie on the elastic line, the relation to be satisfied is 


k'Aej 

~~r~ 



(B-2) 


where 

Nj is the corresponding transition life in the presence of mean 
stress, and the algebraic sign of (positive for tension, negative for 

I 

compression) governs the difference in associated Nj for the two cases. 
Thus, solving for Nj from Eg. (B-2), and substituting into the general 
life relation from [14], 
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where 

Acj is the transition strainrange for basic material under zero 
mean stress, and is given in general terms by Eq. (A-3). 

Nj is the transition life for the basic material under zero mean 
stress, and is given in general terms by Eq. (A-4). 
k is mean stress sensitivity of material, as discussed in 
report, equal to unity according to Morrow's formulation, 
k' is transition life sensitivity of material to mean stress, 
a function of mean stress. 

o^, e^, b and c are the conventional material constants establishing the 
life relationship for the particular material. 

Future experiments are required to establish the k' versus mean 
stress relationship. In the absence of special information, we can 
let k' = 1, which essentially assumes that the cyclic stress-strain 
curve is unaffected by mean stress. 

It is possible to substitue into Eq. (B-4) the expression for 
and Hj from Eqs. (A-3) and (A-4), resulting in 





(B-5) 
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While Eq. (B-5) does not explicitly state how k' varies with it 
is a general framework for future study of mean stress effects. It is 
clear that the life relation applies to both tensile and compressive 
mean stress, depending on the algebraic sign used for a^. 

The life with mean stress can be obtained from Eq. (B-5) 
by following the same inversion procedures outlined in Appendix 
A for the case of zero mean stress. 
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TABLE I. APPLICATION OF THE DOUBLE LINEAR 
DAMAGE RULE TO A COMPLEX MISSION CYCLE 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

hVENT 

LOADING 

CONDITIONS 

NO. CYCLES 
PER MISSION 

CYCLES TO 
FAILURE 

CYCLES TO 
END OF PHASE I 

CYCLES TO 
COMPLETE 

n 

n 


Ae, 

in/in 

^o. 

MPa 

n 

Nf 

Eq.tB-l) 

^^i 

(Eq.2S) 

PHASE II 
Nil 

Eq. (26) 

"i 

"ii 

1 

0.00254 

695 

4 

37180 

23643 

13537 

I.692X lO"'^ 

2.955X lO"'* 

2 

0.00791 

394 

2 

7200 

2474 

4726 

8.084 

4.232 

3 

0.00735 

359 

1 

13650 

6354 

7296 

1.574 

1.371 

4 

0.01017 

268 

6 

5550 

1632 

3918 

36,770 

15.310 

5 

0.00396 

616 

3 

17400 

8873 

8527 

3.381 

3.518 

6 

0.00198 

727 

2 

64000 

45518 

18482 

0.439 

1.082 

7 

0.00848 

172 

1 

33000 

20382 

12618 

0.491 

0.793 

8 

0.01564 

62 

2 

j 2500 

390 

2110 

1 

51.280 

9.479 

9 

0.01045 

3 

1 

' 31325 

19092 

12233 

0.524 

0.817 

10 

0.00932 

66 

1 

42540 

27897 

14643 

0.359 

0.683 

11 

0.01074 

145 

1 

9390 

3706 

5684 

2.698 

1.759 

12 

0.01271 

127 

1 

i 4440 

1122 

3318 

8.913 

1 3.014 

13 

0.01158 

188 

1 

' 4900 

1327 

' 3573 

7.536 

1 2.799 

14 

0.00452 

557 

2 

20605 

11124 

9481 

1 

; 1.798 

2.110 

TOTALS 

I25.539X lo"* 

49.922xI0"'‘ 


MISSIONS TO COMPLETE PHASE I DAMAGE = I/. 0126= 79 
MISSIONS TO COMPLETE PHASE II DAMAGE = 1/.005 = 200 


TOTAL MISSIONS TO FAILURE 


79 + 200 = 279 























TABLE II. APPLICATION OF CUMULATIVE FATIGUE DAMAGE RULES TO THREE STEP 

LEVEL TESTS OF WEBBER AND LEVY (REF. 13) FOR D.T.D. 683 ALUMINUM. 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

(11) 

(12) 

(13) 
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No. 

TEST CONDITION 
% of 9500 cycle block 
at each life level 



CYCLES TO 

FAILURE IN 1000s 


RATIO, PREDICTED/EXPERIMENTAL 

Exp. 

LDR 



Damage 
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Damage^^^ 

Curve 

LDR 

DLDR^*’^ 

n W) 

Damage 

Curve 

mm 


^3 
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36 
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1.32 

1.17 

1.13 
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42 

57 
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801 

672 

668 
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0.92 
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0.75 

3 


IS 

84 

741 

1624 

1151 
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1178 

1292 

2.19 

1.55 

1.59 

4 

4.2 

63 

32.8 

363 

425 

356 

345 

352 

380 

1.17 

0.98 

0.97 

5 

4.2 

42 

53.8 

372 

537 

439 

415 

437 

456 

1.44 

1.18 

1. 17 

6 

4.2 

15 

80.8 

550 

813 
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572 

646 

693 

1.48 

1.16 

1.17 

7 

10 

63 

27 

275 

288 

248 
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238 

257 

1.05 

0.90 

0.86 

1 

8 

10 

42 

48 
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336 

289 

271 
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299 

1.19 

1.03 

1.03 

9 

10 

15 

75 

269 

1 427 

371 

340 

1 

371 

385 

1.59 

! 1.38 

1 

1.38 

Averages 

1.37 

1 1.12 

1. 12 


(a) N^ = 51,900 

N 2 = 414,140 

= 13,800,000 


(b) Using N^ and N^ for calculation of kneepoint 

(c) Using Nj and N^ for calculation of kneepoint 

(d) Using exponent of 0.40 

(e) Using exponent of 0.33 

























Figure 1. - Schematic of damage curve concept of 
summing cumulative damage in complex loading. 
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Figure 2. ~ Cycle ratio relationship for two level 
test as deduced from damage curve, and re- 
placement of damage curves by two straight 

line segments. 
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Figure 3. - Determination of exponent in 
Eq. (14) used in calculating the kneepoint 
for the double linear damage rule. 
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Figure 4 - Determination of coefficient B in 
Eq. (14) used in calculating the kneepoint for 
the double linear damage rule. 
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Figure 5. - Two load level test results for maraging 300 CVM steel showing comparison with predictions by the damage curve approach 
and the double linear damage rule. Data from Ref. (3) and previously unpublished NASA results. 
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Figure 6. - Two load level lest results for SAE 4130 steel showing 
* comparison with predictions by the damage curve approach and 
the double linear damage rule. Data from Ref. (3). 
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Figure 7. - Two load level test results for Ti-6AI-4V with and with- 
out mean stresses showing comparison with predictions by the 
damage curve approach and the double linear damage rule. Data 
from Ref. (12). 
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Figure 8. -Two load level test results for British aluminum 
alloy D.T. D. 683 showing comparison with the damage 
curve analysis and the double linear damage rule. Data 
from Ref. (13). 



Figure 9. - Low cycle fatigue life relations for 
Phase I, Phase II, and total life for Ti-6AI-4V 
at 77° C (170° F). Zero mean stress condition. 
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Figure 10. - Lack of correlation between knee- 
point coordinate and the concept of a discrete 
crack initiation event. 



1. Report No. 2. Government Accession No. 

NASA TM-81517 


4. Title and Subtitle 

PRACTICAL IMPLEMENTATION OF THE DOUBLE LINEAR 
DAMAGE RULE AND DAMAGE CURVE APPROACH FOR 
TREATING CUMULATIVE FATIGUE DAMAGE 


7. Author(s) 

S. S. Manson and G, R. Halford 


9. Performing Organization Name and Address 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 

12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, D.C. 20546 



5. Report Date 

April 1980 

6. Performing Organization Code 

8. Performing Organization Report No. 

E-387 


10. Work Unit No. 


11. Contract or Grant No. 


13. Type of Report and Period Covered 
Technical Memorandum 


14. Sponsoring Agency Code 


15. Supplementary Notes 

S. S. Manson, Case Western Reserve University, Cleveland, Ohio 44106; G. R. Halford, 
NASA Lewis Research Center. 


16. Abstract 

Simple procedures are presented for treating cumulative fatigue damage under complex loading 
history using either the Damage Curve concept or the Double Linear Damage Rule. A single 
equation is provided for use with the Damage Curve approach; each loading event providing a 
fraction of damage until failure is presumed to occur when the damage sum becomes unity. 

For the Double Linear Damage Rule, analytical expressions are provided for determining the 
two phases of life. The procedure involves two steps, each similar to the conventional appli- 
cation of the commonly used Linear Damage Rule. When the sum of cycle ratios based on 
Phase I lives reaches unity, Phase I is presumed complete, and further loadings are summed 
as cycle ratios based on Phase II lives. When the Phase H sum reaches unity, failure is pre- 
sumed to occur. No other physical properties or material constants than those normally used 
in a conventional Linear Damage Rule analysis are required for application of either of the two 
cumulative damage methods described. Illustrations and comparisons of both methods are 
discussed. 


17. Key Words {Suggested by Author(s)) 

Metal fatigue; Crack initiation; Crack 
propagation; Cumulative damage; Linear 
damage; Miner Rule:, Mean stress; 

Life prediction 

18. Distribution Statement 

Unclassified - unlimited 
STAR Category 39 

19. Security Classif. (of this report) 

Unclassified 

20. Security Classif. (of this page) 

Unclassified 

21. No. of Pages 

22. Price* 


For sale by the National Technical Information Service, Springfield, Virginia 221G1 



















National Aeronautics and 
Space Administration 

Washington, D.C. 

20546 


Official Business 

Penalty for Private Use, $300 


SPECIAL FOURTH CLASS MAIL 
BOOK 


Postage and Fees Paid 
National Aeronautics and 
Space Administration 
NASA-451 



! 




POSTMASTER: 


If Undeliverable (Section 158 
Postal Manual) Do Not Return 


/ 


X 


I 



